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The system of 4 differential equations in the external invariant satisfied by the 4 master integrals of the general 
massive 2-loop sunrise self-mass diagram is solved by the Runge-Kutta method in the complex plane. The method 
offers a reliable and robust approach to the direct and precise numerical evaluation of Feynman graph integrals. 


The relevance of the higher order calculations 
for the comparison with nowadays precision mea¬ 
surements in high energy physics is well known 
and comprehensively presented by G. Passarino 
in this conference. 

Therefore can be of some interest the exploita¬ 
tion of an alternative method (but still in the con¬ 
text of the integration by part identities and mas¬ 
ter integrals (MI) 0) to the more common direct 
integration method for the numerical evaluation 
of the MI. 

The method uses directly the differential equa¬ 
tions. Starting from the integral representation 
of the MI, related to a certain Feynman graph, 
by derivation with respect to one of the internal 
masses 0 or one of the external invariants |j) 
and with the repeated use of the integration by 
part identities, a system of independent first or¬ 
der partial differential equations is obtained in a 
number equal to the number of the MI (master 
differential equations). 

Enlarging the number of loops and legs grows 
the number of parameters, MI and equations, but 
does not change or spoil the method. 

* Speaker at AC AT’2002, VIII International Workshop on 
Advanced Computing and Analysis Techniques in Physics 
Research, Moscow, 24-28 June 2002. 


To solve the system of equations (analytically 
or numerically) it is necessary to know the MI for 
a chosen value of the differential parameter. To 
achieve that analytically comes out to be the most 
laborious part of the method and often requires 
some external knowledge, like the assumption of 
regularity of the solution in the value. Moreover, 
if the chosen value is a zero for one of the coeffi¬ 
cients of the derivative of the MI in the equations 
(as is always the case in analytical calculations), 
also the first derivative for that value is neces¬ 
sary to work out the numerical solution for dif¬ 
ferent values of the parameter, but this usually 
comes out to be a simpler task (unless poles in 
the limit of the number of dimensions n going to 
4 are present). 

To test the method we have chosen the simple, 
but not trivial, 2-loop sunrise graph with arbi¬ 
trary masses (J0, shown in Fig.jll 

This graph is one of the topologies of the 2-loop 
self-mass and has 4 MI, the other topologies with 
4 and 5 propagators have one more MI each 00 . 

When the sunrise MI are expanded in (n — 4), 
the coefficients of the poles can be known analyt¬ 
ically for arbitrary values of the external squared 
momentum p 2 , while the finite parts satisfy the 
differential equations. From these equations the 
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Figure 1. The general massive 2-loop sunrise self¬ 
mass diagram. 


analytic expressions for their first order expan¬ 
sion were completed at the special points : 

p 2 = 0 ; p 2 = oo; p 2 = —(mi + m 2 + m 3 ) , the 
threshold; p 2 = -(mi+ m 2 — m 3 ) 2 , the pseudo¬ 
thresholds. 

To obtain numerical results for arbitrary values 
of p 2 , a 4th-order Runge-Kutta formula is imple¬ 
mented in a FORTRAN code, with a solution ad¬ 
vancing path starting from the special points, so 
that also the first term in the expansion is neces¬ 
sary. 

The path followed starts usually from p 2 = 0 
and moves in the lower half complex plane of p 2 to 
avoid proximity to the other special points, which 
can cause loss in precision. For values of p 2 very 
close to a special point, we start from the analyt¬ 
ical expansion at that special point. Remarkable 
self-consistency checks are provided by choosing 
different paths or different starting points to cal¬ 
culate the same value. 

The execution of the program is rather fast and 
precise: with an Intel Pentium III of 1 GHz we 
get values with 7 digits requiring times ranging 
from a fraction of a second to 10 seconds of CPU, 
and with 11 digits from few tens of seconds to one 
hour. 

If A = L/N is the length of one step, L is the 
length of the whole path and N the total num¬ 
ber of steps, the 4th-order Runge-Kutta formula 
discards terms of order A 5 , so the whole error is 
= NA 5 = L 5 /N 4 , and a proper choice of L 
and N allows the control of the precision. 

Indeed we estimate the relative error, as usual, 
by comparing a value obtained with N steps 


with the one obtained with N/10 steps, crk = 
[V(./V) — V(N/10)]/V(N), to which we add a cu¬ 
mulative rounding error e cre = \/N x 10 —15 , due 
to our 15 digits double precision. 

Comparisons are done in |q] w ith some val¬ 
ues present in the literature l] with excellent 
agreement. 
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